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Abstract 

Many situations in physics, biology, and engineering consist of the transport of some physical quantity through a network 
of narrow channels. The ability of a network to transport such a quantity in every direction can be described by the average 
conductivity associated with. When the flow through each channel is conserved and derives from a potential function, we 
show that there exist an upper bound of the average conductivity, and explicitly give the expression for this upper bound as 
a function of the channel permeability and channel length distributions. Moreover, we express the necessary and sufficient 
conditions on the network structure to maximize the average conductivity. These conditions are found to be independent of 
the connectivity of the vertices. 



Examples of transport phenomena through a network 
of channels abound in nature and engineering: blood 
flow through the microcirculatory system, water trans- 
port through the venation of a leaf, water and electricity 
supply in a city, heat conduction through an open cell 
material, etc. Therefore, the search for a network struc- 
ture optimizing the transport processes may be useful 
for a better understanding of the structure of natural 
networks and for the conception of optimized materials. 
Different functions can be optimized, such as dissipated 
power, volume or surface area of the channels. Different 
constraints can be imposed: topology, flow rate, volume, 
etc. Various models have been proposed to understand 
the structure of natural networks 0,0, mostly based on 
the assumption of local optimization. The idea is to move 
the position of a given junction, with the other junctions 
and the topology fixed, in order to optimize one of the 
above functions. In a previous work, we considered the 
global optimization of the whole network structure for 
its electrical property . We showed the existence of an 
upper bound for the average electrical conductivity of a 
network made of uniform wires, for a given amount of 
material per unit volume, and proved two necessary and 
sufficient conditions on the structure of the network for 
the upper bound to be reached. We then used these re- 
sults to derive the expression of electrical conductivity of 
dry foams. In the present paper, we generalize the idea of 
an optimal structure for transport processes to networks 
made of non-uniform pipes and to any flow process, and 
derive three more general necessary and sufficient condi- 
tions for maximizing the transport properties. Surpris- 
ingly, we find that these conditions do not depend on the 
connectivity of the junctions. 

Let us first denote each pipe by a pair of indices (i,j) 
corresponding to the labels of its two ends. We con- 
sider pipes for which the aspect ratio is very large, so 



we can univocally define for each pipe a length lij, 
a local cross-sectional area s«(/), and a local permeabil- 
ity kij (I) (kij (I) is function of Sy (I) and both can vary 
with the curvilinear coordinate I along the channel) . Let 
us assume that the flow through each channel is directly 
related to the gradient of a potential function V, so the 
flow vector I, the permeability and V are related by: 
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(i) 



Moreover we suppose that the flow is in a steady state (so 
V satisfies Laplace equation: V 2 V = 0), and is conserved 
through each pipe and each junction. We then define 
the " dissipated power''' associated with the flow in the 
network as: 



I • Wdl. 



(2) 



Using the assumption of steady state flow, we can rewrite 
the dissipated power as: 



P — Uijiij, 

(id) 



(3) 
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where Uij and iij are respectively the potential differ- 
ence and the flow rate through the channel Since 
the flow is conserved through each pipe and each junc- 
tion, the total dissipated power in the network is equal 
to the product U.I, where U is the potential difference 
between the inlet (source) and the outlet (sink) of the 
network and I is the crossing flow rate, as for the dissi- 
pated power of an electrical network. Pursuing the anal- 
ogy, it follows that the actual distribution of flow rates 
for a given total flow rate is the one which minimize the 
dissipated power, and subsequently the monotonicity law 
of Rayleigh is valid Q. 0. 0| : If any of the resistances of 
a circuit are increased, the effective resistance between 
any two points can only increase. If they are decreased, 
the effective resistance can only decrease. Although the 
monotonicity law cannot predict the sign of the effective 
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resistance variation when some resistances are increased 
and some others are decreased (the variation obviously 
depends on the structure of the network) , it will be help- 
ful to determine the optimal geometry of a network for 
its transport properties. 

A. Maximal average conductivity of a network 

On a macroscopic scale, the network can be seen as 
a continuous medium, a priori anisotropic, whose trans- 
port properties are described by an effective conductivity 
tensor a - the conductivity being defined as the perme- 
ability per unit length (2D network) or per unit area (3D 
network). The average conductivity associated with this 
tensor is defined as: 




outlet 



FIG. 1: Schematic representation of a network shorted with 
parallel sheets of infinite conductivity. The resistance of the 
truncated channel is equal to Ax/ (fcij (I) cos Oiij). 



a m = (u.g.u) , (4) 

where the brackets indicate that the term inside is av- 
eraged in all the directions of the unit vector u. This 
parameter can be simply related to the trace of a and 
the dimension of space d by: 

a m = ^Tr [g] . (5) 

The average conductivity characterizes the ability of the 
network to transport the physical quantity associated 
with in all the directions. We shall show the existence 
of an upper bound for <j m , which can be expressed as a 
function of the channel permeability and channel length 
distributions: 

L d is the hypervolume of the network on which the con- 
ductivity tensor is defined, and k%j is the average perme- 
ability of the channel defined by: 

kij = -r [ %3 kij^dl. (7) 

Hj Jo 

We present here the demonstration for a three- 
dimensional network (d = 3). We first study the case 
of a network made of straight pipes: imagine that this 
network is shorted with thin parallel sheets of infinite 
conductivity, perpendicular to the direction x of the ap- 
plied potential difference, and separated from each other 
by Ax, as in FigQ] (Ax has to be small compared to the 
typical channel length, but large compared to the typical 
channel diameter, so the flow lines in a given channel are 
mostly along its axis). According to the monotonicity 
law, the conductivity of this shorted network is higher 
than the conductivity of the original one. Furthermore 
the potential is uniform on each sheet and the resistance 
AR{x) of the network slice at position x corresponds to 
the parallel association of the truncated channels that it 



contains. Since the resistance corresponding to a given 
truncated channel is equal to Ax/ (fey (I) cosoiij), 

where cty is the angle between the channel (i,j) and the 
axis x (the angle is defined such that cos ay > 0), we 
have: 

1 - V P(x-x- x ^ MQcos"ij (8) 
AR(x) ^ [,t,j) Ax ' [) 
(hi) 

where the sum is carried out on all the channels of the 
network by introducing the function P(x; Xi, Xj) which 
take the value 1 if the channel is intersected by the 
equipotential plane passing by x (i.e. if x is between Xi 
and xj) and otherwise. The total resistance is given 
by the sum of the slice resistances. Using the fact that 
the product of the average of a set of positive values 
/2j •••> /at} by the average of their inverses is always 
greater or equal to one: 

and passing to the continuum limit, we obtain a lower 
bound for the resistance of the shorted network: 

^ /*oo dx 
— < I y^P{x,Xi,x j )ki j {l)cosa ij - F1 . (10) 

We can switch the sum and the integral of this expression, 
and it follows after integration that: 

°i S) ^ 73 £ ^ ^ cos2 a ^ ' ( 1 1 ) 

where the conductivity of the shorted network in the x Ax- 
is) j 

rection is defined by a x — L ^ r > an d L z being the 
network lengths in directions y and z, and L 3 — L x L y L z 
the volume on which the conductivity is defined. Using 
the same arguments in the two other directions and the 
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fact that the sum of the squared direct cosines is equal to 
one, we see that inequality is true for a network made 
of straight channels. It is clear that the average con- 
ductivity of a network with curved channels is bounded 
too: in this case, we can built a new network by keep- 
ing the junctions fixed in their positions but linked with 
straight channels. Let & be the length of the straight 
channel Moreover we can choose its local perme- 

ability k'ijil) to be equal to the local permeability of the 
corresponding curved channel on an arbitrary portion l'^ 

So 



of its length, say: fey (I) = feij(Z) for < / < 



di 



the resistance of the straight channel, equal to J Q - y ^ , 

is lower than the original curved channel. We know the 
conductivity a' m of such a network is bounded: 



< 
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(12) 



On the one hand the conductivity of this network is 
higher than the conductivity a m of the original network 
(from the monotonicity law), and on the other hand 
I'ijk'ij < hjkij, so inequality JSJ) holds for a network made 
of curved channels as well. 



B. Optimizing transport 

In addition to the existence and the expression of an 
upper bound for the average conductivity, we show that 
the average conductivity reaches the upper bound if and 
only if the following three conditions are satisfied: 

a) Each channel has a uniform local permeability 
along the channel: kij(l) 



h ■ — h ■ 



b) All the channels are straight. 



c) Every junction (i) between channels satisfies 
J2j kij e ij — 0, where ey are outward-pointing unit 
vectors in the directions of adjoining channels. 



Note that the last condition is equivalent to a force 
balance equation at each vertex, the weight of the force 
pulling along a channel being analogous to its permeabil- 
ity. Furthermore, it is worth noting that the three condi- 
tions are independent of the connectivity of the junctions. 
As an illustrative example of this property, the periodic 
square, hexagonal and triangular networks built with a 
same set of channels have the same conductivity. 

Necessity of the conditions: Let us suppose that 
the conductivity of the network is equal to its maximal 
value: 



d ^ 



(13) 



and let us make some infinitesimal changes in the network 
structure. First, suppose we vary the local permeability 



of a given channel by a small amount 6kij(l), with- 
out altering its length or its average permeability fey , 
so: 



6k i3 {l)dl = 0. 



(14) 



Since Zy and kij remain constant, the corresponding vari- 
ation 5a m of the conductivity is zero. But obviously we 
can choose a variation Skij (I) such that the resistance of 
the channel is increased: 



(15) 



In order for the monotonicity law to be satisfied, the 
variation of the resistance has to be zero, which can be 
mathematically expressed as: 
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+ A 6k i:j (l)dl = 0, 



(16) 



where A is a Lagrange multiplier. Since this equality has 
to be true for any variation 8kij(l), it follows that kij (I) 
has to be uniform along the channel, and so condition a. 
is indeed required. 

Necessity of conditions b) and c) is proved by means of 
arguments similar to those used in our previous workQ: 
first, imagine that we change the length of a given channel 
of a network for which condition a) is satisfied, the 
positions of all junctions and the lengths of all the other 
channels remaining unaltered. To this variation of length 
SUj corresponds a variation of the conductivity: 



6a r , 



^ d kijSlij, 



(17) 



implying that Sa m and 5Uj have the same sign. However, 
if the length of the channel is increased (Shj > 0), the 
resistance of the channel is increased too, and it follows 
from the monotonicity law that the conductivity can only 
decrease (5<r m < 0). So the variation of the length has 
be zero at first order, implying the necessity of condition 

Now, imagine a network for which conditions a) and b) 
are satisfied. In the previous section, we showed the exis- 
tence of an upper bound of the conductivity in the direc- 
tion x for a network made of straight channels, by using 
two successive inequalities: first, the conductivity of such 
a network is lower than the conductivity of the same net- 
work intersected with zero resistance sheets; second, the 
shorted network conductivity itself is bounded, using the 
fact that the equivalent resistance of iV resistive elements 
in series arrangement is greater or equal to the equiva- 
lent resistance of the same resistive elements in parallel 
arrangement times iV 2 (inequality ©). So in order to 
get the exact upper bound, these two inequalities have 
to become strict equalities. The first one implies that 
the presence of sheets does not modify the distributions 
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of potentials, and so the potential in the channels is a 
function of x only. To see this, increase progressively the 
resistance of the sheets up to infinity (what corresponds 
to the initial network). From the monotonicity law, this 
can only decrease the conductivity of the network. The 
only way for the conductivity to stay at its maximum 
value when increasing the resistance of each sheet is by 
having no current through them, and so the potential in 
the channels is a function of x only. The second equal- 
ity requires that the resistance of every slice (of equal 
thickness) is the same, or equivalcntly, the resistance of 
a slice of arbitrary thickness x is simply proportional to 
x. Hence the potential is indeed linear in x. Examina- 
tion of flow conservation at a vertex in such a potential 
immediately leads to condition : 

kjj cos ajj.sgn (xj — Xj) = 0, (18) 

3 

(the term sgn [xi — Xj) is introduced in order to satisfy 
cosay > 0). This equality is nothing but the projection 
of condition c) on the axis x. Since the same argument 
holds in the two other directions, we prove that condition 
c) is indeed required. 




FIG. 2: Netted venation of a leaf. The angles between adja- 
cent veins are correlated to their cross-section areas Hell. 



Sufficiency of the conditions: Now consider a net- 
work for which conditions a), b), c) are satisfied, and 
suppose that a potential difference U x is applied between 
the two regarding faces orthogonal to x. Let us show 
that the trial potential function defined as <p = —jf-x 
is the physical solution. We first check that the flow is 
conserved at each junction under the stated conditions: 
the current in the straight channel is given by: 

J-'x 

so: 

J y = ~r ex ■ H / '-.' <> ',' = °- ( 2 °) 



The trial potential function also satisfies the boundary 
conditions, and so is the correct physical solution. Now 
we can show that the average conductivity is equal to 
the upper bound. The potential is uniform on planes 
perpendicular to the direction x, and so the system is 
unaltered when intersected by thin parallel sheets of in- 
finite conductivity orthogonal to that direction. We pre- 
viously calculated the elementary resistance AR(x) of a 
slice of thickness Ace at position x for such a network 
(Eq. [SJ. The global flow rate I x , the elementary resis- 
tance AR(x), and the elementary potential difference A</> 
across the slice are related by: 

As Ax 

But here both the potential gradient and the global cur- 
rent are independent of x, and so is -^7. Using Eq. (jHJ) 
and integrating in the x direction, we finally get the ex- 
pression of the conductivity in the x direction: 

° x = Jd X! cos2 ai i- ( 22 ) 

L m 

This expression corresponds indeed to the upper bound of 
the conductivity along x, as expressed in The same 
reasoning can be applied for the conductivity in the two 
others directions, and so the sufficiency of conditions a), 
b) and c) is proven. 

Let us apply the preceding results to two very com- 
mon flow profiles, namely plug flow and Poiseuille flow. 
Many transport phenomena through pipes are described 
by plug flow, such as the flow of fluids in porous con- 
ducts, the flow of heat by conduction, the electrical cur- 
rent, or diffusive flow. All these flows are engendered 
by a gradient of a potential function (such as pressure, 
temperature, electric potential, material concentration), 
and the permeability for all these situations is directly 
proportional to the cross-sectional area. Moreover, if 
the conductivity has the same value do for all pipes (so 
hj = o'oSij), the upper bound of the average conductivity 
simplifies to 

0>n < "Woe, (23) 
a 

where e is the volume fraction of the continuous phase. 
As a consequence, the maximal conductivity does not 
depend on the permeability and length channel distribu- 
tions, but only on the total amount of material that the 
network contains (there exists a universal upper bound 
of the quantity (cr m /e), the average conductivity per unit 
volume of material) . Note that the expression of the up- 
per bound is equivalent to the Hashin-Shtrikman bound 
for the electrical conductivity of an heterogenous mate- 
rial in the limit of small volume fraction of the conducting 
phase 0||]. 

Condition c) leads to different optimal structures de- 
pending on the kind of flow profile in the channels. In 
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the case of a plug flow profile with idential conductivity 
(To, this condition becomes: 

^8^3=0. (24) 
j 

In the case of Poiseuille flow, the permeability varies like 
the square of the cross-sectional area. Thus, if we assume 
again the same conductivity <jq for every pipe, condition 
c) becomes: 

E4 e ^ = °- ( 25 ) 

j 

As a concluding remark, it is worth noting that some 
natural networks, like leaf venation[Iol ITT| . have a very 
well defined structure, in the sense that the angles 



between adjacent veins are correlated to their cross- 
sectional areas (see Fig. [2J. This fact presumably cor- 
responds to some optimization principle. The conditions 
for transport optimization presented here may give an 
explanation for these typical patterns. In the particu- 
lar case of leaf venation, veins are composed of a com- 
plex tangle of smaller tubes , hence we may expect the 
relation (|24|l rather than l|25|> to be satisfied (although 
veins are not fully impermeable) . Alternatively, this cor- 
relation between angles and cross-section areas could be 
explained by the o ptim ization of the mechanical stabil- 
ity of the leaf 0, Il3| . Experimental study of the leaf 
venation structure is subject to a current investigation 
to be compared with the transport optimization criteria 
presented in this paper. 

M.D. thanks H. A. Stone for encouraging this research. 



[1] M. Zamir, J. theor. Biol. 62, 227-251 (1976). 

[2] M. A. Changizi and C. Cherniak, Can. J. Physiol. Phar- 
macol. 78, 603-611 (2000). 

[3] M. Durand, J.-F. Sadoc, and D. Weaire, Proc. R. Soc. 
Lond. A. 460, 1269-1285 (2004). 

[4] J. W. S. Rayleigh, On the theory of resonance. In Col- 
lected scientific papers, volume 1, 33-75 (1899). 

[5] J. H. Jeans, The mathematical theory of electricity and 
magnetism, 5th edition, Cambridge University Press 
(1925). 

[6] A priori cr m can be a function of the trace and the deter- 
minant of <t, which are the two invariants of a symmetri- 
cal 2 nd order tensor, but using the additivity property of 
tensors and dimensional analysis, it follows that a m can 
only depend on the trace, in a linearly way. The prefactor 



1/d comes naturally by studying the case of an isotropic 
network. 

[7] Z. Hashin and S. Shtrikman, J. of Appl. Phys. 35, 3125- 
3131 (1962). 

[8] S. Torquato, Random Heterogeneous Materials, 

Springer- Verlag, New- York (2002). 
[9] A. H. de Boer and V. Volkov, Plant, Cell and Environ- 
ment 26, 87-101 (2003). 

[10] S. Bohn, B. Andreotti, S. Douady, J. Munzinger, and Y. 
Couder, Phys. Rev. E 65, 061914 (2002). 

[11] Y. Couder, L. Pauchard, C. Allain, M. Adda-Beddia, and 
S. Douady, Eur. Phys. J. B 28, 135-138 (2002). 

[12] A. Roth-Nebelsick, D. Uhl, V. Mosbrugger, and H. Kerp, 
Annals of Botany 87, 553-566 (2001). 

[13] M. Durand (unpublished). 



5 



